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ABSTRACT 

We re-estimate the surface magnetic fields of neutron stars (NSs) in Be X-ray bi¬ 
naries (BeXBs) with different models of torque, improved beyond Klus et al. (2014). 
In particular, a new torque model (Dai &: Li 2006) is applied to three models of mag¬ 
netosphere radius. Unlike the previous models, the new torque model does not lead 
to divergent results for any fastness parameter. The inferred surface magnetic fields of 
these NSs for the two compressed magnetosphere models are much higher than that for 
the uncompressed magnetosphere model. The new torque model using the compressed- 
magnetosphere radius (Shi, Zhang & Li 2014) leads to unique solutions near spin equi¬ 
librium in all cases, unlike other models that usually give two branches of solutions. 
Although our conclusions are still affected by the simplistic assumptions about the 
magnetosphere radius calculations, we show several groups of possible surface magnetic 
field values with our new models when the interaction between the magnetosphere and 
the infalling accretion plasma is considered. The estimated surface magnetic fields for 
NSs BeXBs in the Large Magellanic Cloud, the Small Magellanic Cloud and the Milk 
Way are between the quantum critical field and the maximum “virial” value by the spin 
equilibrium condition. 
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— X-rays: binaries 


^National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 

^College of Material Science and Chemical Engineering, Hainan University, Hainan 570228, China 

®Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 
210093, China 

^Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
100049, China; zhangsn@ihep.ac.cn 

^Department of Astronomy, Nanjing University, Nanjing 210093, China 



- 2 - 


1. Introduction 

Be X-ray binaries (BeXBs) comprise a neutron star (NS) and a fast rotating Be-type main 
sequence star and are the most numerous class of high-mass X-ray binaries (HMXBs; Liu et ah 
2006). Generally, a 8Mq — 18Mq Be star in a BeXB is surrounded by a circumstellar disk (Knigge 
et al. 2011; Reig 2011), where Mq is the mass of the sun. A circling NS accretes matter from the 
disk arising from the wind of the Be star, when it passes through the circumstellar disk and so 
most BeXBs are transient systems due to the eccentric orbits (e > 0.3) (Reig 2011). 

The accretion process of an NS in a BeXB has been discussed by many researchers such as 
Shakura et al. (2012) and Ho et al. (2014). For X-ray binaries (XBs), a necessary condition for the 
formation of an accretion disk is that the specific angular momentum of gravitationally captured 
matter is sufficiently high, or the accretion wind will persist (Burnard et al. 1983; Shakura et al. 
2012, 2014). Reig (2011) discussed that Be stars as fast rotators might eject photospheric matter 
with sufficient energy and angular momentum due to some weak processes, such as gas pressure or 
non-radial pulsations, so that a Keplerian disk is formed. Porter & Rivinius (2003) considered that 
emission lines and infrared excess in BeXBs might originate from an equatorial disk. James (2010) 
listed 14 HMXBs with quasi-periodic oscillations, and many researchers considered that QPOs in 
XBs were an indication of an inner accretion disk around an NS (Miller et al. 1998; Stella & Vietri 
1998; Osherovich & Titarchuk 1999; Abramowicz & Kluzniak 2001; Zhang 2004; Li & Zhang 2005; 
Shi & Li 2009, 2010; Shi 2011; Shi, Zhang & Li 2014, etc.). 

Giacconi et al. (1973 ) first observed the pulse period’s variation of Hercules X-1 and later 
two states (spin-up or spin-down) of NSs were observed (Nagase 1989; Bildsten et al. 1997). Many 
authors (such as Ghosh & Lamb 1979; Wang 1995; Kluzniak & Rappaport 2007) proposed different 
physical mechanisms in order to explain the change of the spin period of the NS in an XB. The 
changing spin period (P) of the NS in an XB was attributed to the change of the angular momentum 
of the NS due to accretion or outflow. 

In most views, the corotation radius is considered as the dividing line between the spin-up 
state and the spin-down state for NSs in XBs. In the early literature such as Ghosh & Lamb 
(1979), the accretion matter from the donor star rotates around a central NS utill it reaches the 
inner disk terminal radius (n), where the last angular velocity is the same as the angular velocity 
of the NS. The accretion disk has always been considered to be truncated at the magnetosphere 
radius (rm), where the pressure of the accreting material is compensated by the magnetic pressure 
of the magnetosphere (Pringle & Rees 1972). Ghosh & Lamb (1979) supposed that the magnetic 
field lines were threaded in the Keplerian accretion disk in a broad transition zone, which has been 
extended to derive the torque acting on the NS. 

Illiaronov & Sunyaev (1975) proposed that accretion could be continued only when the ac¬ 
creting pulsar’s magnetosphere rotates more slowly than the Keplerian angular frequency (12^ = 
y^GM/r^) of the accretion matter at rm, where G, M, r are the gravitational constant, the mass of 
an NS, the radius to the center of the NS, respectively. In other word, if the magnetospheric radius 
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lies outside the corotation radius (vco = ^/GMP^/47r^), at which the Keplerian angular frequency 
of the accretion matter is equal to the NS spin angular frequency (fig), the accretion would be 
inhibited, i.e., the propeller effect emerges in the system (Shvartsman 1970). The remaining state 
of the NS (neither spin-up nor spin-down) is called spin equilibrium and it was believed to emerge 
when Tco = Tm (Pringle &: Rees 1972, Klus et al. 2014). 

Klus et al. (2014) obtained the spin period, change rate of spin period (P), and X-ray lumi¬ 
nosity for 42 BeXBs with 2 states (spin-up or spin-down) in the Small Magellanic Cloud (SMC). 
An unexpected conclusion was obtained: the magnetic fields of the NSs in BeXBs are either much 
stronger or lower than those in low mass X-ray binaries (LMXBs). Ho et al. (2014) continued 
to discuss the two different kinds of solutions and they concluded that the strong magnetic field 
solutions from the non-spin equilibrium condition are more compelling because these solutions are 
close to the solutions from the spin equilibrium condition. They believed that spin equilibrium of 
BeXBs finally united different neutron star populations. 

This paper is organized as follows. In Section 2, we review five available torque models and 
present a new torque model. In Section 3, we estimate the surface magnetic fields of NSs in BeXBs 
with these torque models and by assuming spin equilibrium for disk accretion. In Section 4, we 
give a discussion on the surface magnetic field of NSs in BeXBs. In the last section, we summarize 
this research. 


2. Models on dimensionless torque 
2.1. Five available models of dimensionless torque 

The changing spin period of an accretion X-ray pulsar was discovered for a long time ago and 
it is widely accepted that the accretion torque acting on a NS system leads to the change of the 
spin period, i.e., 

-P = iVPV2vr/efr , (1) 

where N is the total torque acting on an NS system and Jeff is the effective moment of inertia 
of the NS (Ghosh et al. 1977). Ghosh &: Lamb (1979) first made a magnetically threaded disk 
model (hereafter GL) to explain the changing spin period of NSs in XBs. Wang (1995) obtained 
three kinds of dimensionless torques in different conditions. Kluzniak &: Rappaport (2007) obtained 
the results on torques acting on the central star and the disk luminosity produced by the released 
gravitational potential energy and the mechanical energy input to the disk via the magnetic torques; 
those results are independent of the viscosity prescription. We will review the main results of the 
models below. 
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2.1.1. The Model of Ghosh and Lamb 


Ghosh & Lamb (1979) considered a magnetic NS with the dipolar magnetic field accreting 
matter from a Keplerian disk. The steady accretion flow has an axial symmetry configuration and 
the magnetic moment is aligned with the rotation axis of the NS. They found that the magnetic 
coupling between the NS and the plasma outside ri was extremely important for the physical 
mechanism of the changing spin of the NS in accretion magnetic XB. The inner disk terminal 
radius (n) in GL can be written as follows 


n Ri 0.52rA 


0.52(- 

2GMM2 



( 2 ) 


where ta is the Alfv& radius for spherical accretion, M the accretion rate, and the magnetic 
moment of an NS. 


They expressed the total torque N = Nq + Nmag, where Nq is the material torque acting on 
the NS system when all the angular momentum of the accretion matter is transferred to the NS 
and Nmag is the magnetic torque produced by the magnetic coupling. The two variables can be 


expressed as follows 


and 


No = M^/GMn , 


-^mag — 


r^out 

/ r^B^B^dr , 
J r\ 


(3) 


(4) 


where rout is the outer radius of the transition zone in the model of Ghosh & Lamb and it is always 
substituted by inhnity because rout The toroidal magnetic held {B^p) is supposed to be 

generated from the vertical magnetic held {B^) due to the differential rotation between an NS and 
its circling disk. The vertical magnetic held is supposed as a dipole magnetic held (i.e. B^ = /t/r^) 
in all the models except Wang’s models below, where jj, is the magnetic moment of an NS. 


The total torque can also be expressed as follows 


N = n * No , 

where n is a dimensionless accretion torque and it can be expressed approximately by 

n = 1.39{1 - w[4.03(l - - 0.878]}(1 - , 

where oj is the fastness parameter dehned as 

= (-N. 




OJ = 


^kin) 


(5) 

( 6 ) 

(7) 


Equation (6) is accurate to within 5% for 0 < w < 0.9 and the fastness parameter can be in the 
range 0 < a; < 1 in GL if the inner radius could be extended all the way to the center of the NS, 
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but the real inner radius must be more than the radius of an NS and this restriction is adopted in 
this study below. 

There are several problems in the model. (1) Wang (1987, 1995) pointed out that in GL the 
Btp component was overestimated so that the magnetic pressure exceeds the gas pressure in the 
outer region of the disk. (2) The zero torque will emerge when oj ~ 0.349 or a; ~ 0.987, which 
means that the torque reverses twice between a positive torque and a negative torque when the 
accretion rate increases from zero to a certain value (corresponding to cu < 0.349), i.e., two spin-up 
and spin-down cycles should be observed during a monotonic change of luminosity over a large 
range. However, such a phenomenon has not been observed so far. 



Fig. 1.— Relations between the dimensionless torque and the fastness parameters in five models. 
The dotted lines show the fastness parameters corresponding to zero torque. 


2.1.2. Three models of Wang (1995) 

Wang (1995) gives a relation between the azimuthal magnetic field and the generated toroidal 
magnetic field due to the differential rotation as follows 

^ = 7r<^(ns - f^k) , (8) 

where 7 is a numerical parameter and is the dissipation timescale of B^. Then Wang (1995) 
discussed three conditions that the dissipation timescale is determined by the reconnection outside 
the disk, the turbulent diffusion in the disk, and the Alfven velocity, respectively. The three models 
in the three conditions are called Wang-(A), Wang-(B), Wang-(C) sequentially below. The relations 
between Bp and B^, can be simplihed when all the characteristic parameters are adopted (Wang 
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1987, 1996) as follows 

(A) ; B^p/B^ = ris/llk — 1 for fig < ^kj 1 ~ for fig > 

, (B) : Bp/B^ = flg/flk - 1 , (9) 

_ (C) : \Bp/B,\ = , 

where p is the thermal pressure of the gas in a standard a disk (Shakura &: Sunyaev 1973). The 
thermal pressure can be expressed as 

p = 1.23 * - 1.12r(rfo2)4/5gpj^-ig-2 ^ 

where R denotes the radius of an NS, L the X-ray luminosity; and the subscripts 6, 37 and 1.40 
express the quantities in units of 10® cm, 10^^ ergs“^, and 1.4 times the mass of the Sun, respectively. 
A negative vertical magnetic field (i.e. B^ = —rjp/r^) is adopted to correct the direction of Bp in 
these models, where r/ is a screening coeffeicient as a constant {rj < 1, Ghosh & Lamb, 1979) and 
the characteristic value r/ = 1 is adopted in this study. The surface magnetic field of the NS in 
a BeXB would be amplified by a factor of about 2 when the extremely small value rj = 0.2 was 
adopted. The inner boundary of the Keplerian disk in Wang’s models (inside Tco) can be obtained 
according to 

{BpB,)\r=r, = -MVGMn/{2rf) . (11) 

The fastness parameter is restricted as 0 < cu < 1 in the above Equations (9, 11) and we can 
simplify Equation (9) at the inner radius of Keplerian disk as follows 

(A) : cj - 1 , 

(B) : uj - 1^_ (12) 

(C) : \r=ri ■ 




r=r-. — \ 


The dimensionless torques can be written as follows 


(7/6)-(4/3)cj-Kl/9)i.j^ 

1—CJ 

(7/6)-(4/3)c,; 

1—CJ ’ 


(A): 
n = < (B) ; 

. 1 I 1 (.j 57/40 .H(l_y)l/2 POO (j/-l)l/2 

( • J- + 3 (1-01/2 Uo; 97/40 «y-Jl 97/40 “2/1 ' 


(13) 


2.1.3. The Model of Kluzniak and Rappoport 

Kluzniak & Rappaport (2007) also discussed the properties of an accretion disk surrounding a 
magnetized star and they first incorporated another component of the material torque (— 
omitted in other models, where R is the radius of an NS (hereafter we will mark the model with 
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KR). That component can be considered as a torque from the remaining angular momentum of the 
accreted plasma after the plasma falls on the surface of the NS and corotates with it. 


This model predicted a dimensionless torque for = Bz{l — ^k/^s) inside Tco as follows 


^ _ 10 


1 u 


9 ISl-cj 

The relation between oj and r\ can be obtained as follows 




(14) 


(15) 


where the fastness parameter must be less than unity. 


As shown in Figure 1, the relation between the dimensionless torque and the fastness parameter 
is very similar to Wang (A) and Wang (B). KR also has a very similar fastness parameter for the 
zero torque with Wang (C). 

All of the above models imply that the fastness parameter is less than unity and it means that 
the accretion process is not in a propeller state. However, in reality the fastness parameter can 
also be slightly higher than unity because the accretion in the propeller phases has been observed 
(Finger et al. 1994). As shown in Figure 1, the dimensionless torques in the above models are 
divergent when cu —)• 1 (also see Equation (29) in Wang 1987), i.e., the models result in an infinite 
torque. The reason is that, e.g., Wang (1987) assumed that within the boundary layer of the 
accretion disk the toroidal field changes with the same law as that on the accretion disk, which is 
not verified. 


2.2. A New torque model 

Besides the torque equilibrium condition (i.e. the method of balancing the magnetic and 
viscous torques), we can also estimate the inner disk terminal radius by the pressure (or energy) 
equilibrium condition. Both theoretical and numerical investigations (e.g., Lynden-Bell & Boily 
1994; Lovelace et al. 1995; Uzdensky et al. 2002) have shown that when the toroidal field becomes 
comparable to the poloidal field, the fields open up and the star-disk linkage is broken. Therefore, 
in the boundary layer the magnetic pitch is always around unity. Inserting this condition into 
Equation (28) in Wang (1987), one can derive the expression of the magnetospheric radius, which 
is the same as that in the case of pressure/energy balance. This is exactly what Wang did in his later 
paper (Wang 1996) to estimate rm by torque balance and the divergence problem can be avoided 
naturally. Numerical calculations by Long et al. (2005) also confirmed that the magnetospheric 
radii in both cases coincide with each other under this condition. Dai & Li (2006) thus adopted 
the inner disk terminal radius by the energy equilibrium condition and made a new torque model. 
In this subsection, we extend this new model to the inner disk terminal radius considered as three 
kinds of magnetosphere radii from the pressure (or energy) equilibrium condition. 



2.2.1. Basic prescriptions 


We adopt the typical assumptions of NSs in BeXBs as follows. (1) The spin of an NS with a 
dipolar magnetic field is aligned to the magnetic axis and perpendicular to the Keplerian accretion 
disk. (2) The azimuthal component of the magnetic field is generated by rotation shear when the 
poloidal component keeps the dipolar field unchanged. (3) The magnetosphere is force-free and 
reconnection takes place outside the disk (Wang 1995). (4) The interaction between the NS and 
the accretion disk mainly comes from the magnetic field anchored on the NS surface. (5) The 
accretion process is steady. (6) The accreted matter at the magnetosphere radius corotates with 
the NS and the magnetosphere, i.e, the accretion matter keeps the last angular momentum at the 
magnetosphere radius. Here the magnetosphere radius (rm) is considered to be the inner radius of 
the Keplerian dsik. The total material torque acting on the NS system can be described as 

Nms,t = M\/GMr^-Mnsrl^ = M^/gWt^{1-uj) , (16) 

and so the total torque can be expressed as follows 


1 ^ — Nmat + Njnag 


(17) 


Wang (1995) found 

^ = I - 1), Hs < Hk , 

B, \ 7(1 - Hk/a), a > Hk , ^ ^ 

where 7 is a constant and hereafter we suppose 7 = 1 , which means the highest steepness of the 
transition between Keplerian motion inside the disk and the corotation with the NS outside the 
disk (Wang 1995; Dai & Li 2006). The magnetic torque in Equation (4) can be integrated by 
combining Equation (17) and we obtain 


Nmag — 


-2co + cj^), ca < 1 , 

- 1 ), cj > 1 , 


(19) 


where fi = BE? is the magnetic moment and the original expression in Dai & Li (2006) in the 
parentheses of the above first sub-equation was 1 — 2a; -|- |a;^ due to typo. 

Combing Equations (5), (16), (17), and (19), we obtain 


n = 


., 2^-3 


(1 w) -I- 3 

(l_^) + i^-^ 


3 My/GMri, 


(^ — 2a; o;^), a; ^ 1 , 

(|a;“^ — 1), a; > 1 . 


( 20 ) 


2.2.2. Three Models of Magnetosphere Radius 

Lamb et al. (1973) discussed the spherical accretion process onto a compact star and gave 
a definition of the magnetosphere radius. Many authors also explored the outer boundary of the 
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magnetospheres of pulsars (e.g., Eisner & Lamb 1977; Burnard et al. 1983; Mitra 1992; Li & Wang 
1995; Weng & Zhang 2011) for the disk accretion process and they obtained analogous conclusions 
according to the spherical accretion. The general expression of that radius in those models can be 
written as follows 

rm = krA ■ (21) 

Klus et al. (2014) adopted Alfv& radius, which can be expressed as follows 

rmi = ta = 3.69 * cm , (22) 

where B denotes the surface magnetic field of an NS; and the subscripts 12 express the quantities 
in units of 10^^ G . 


Shi et al. (2014) extended the definition of the magnetosphere radius and they considered that 
the magnetohydradynamic equation, including gravity, barometric pressure, magnetic pressure, and 
inertial force acting on the accretion plasma at the magnetosphere radius in a corotation frame of 
reference, must be satisfied. They considered that the magnetic field is compressed by the accretion 
matter utill the accretion matter reaches the magnetosphere radius but the dipolar magnetic field 
inside the magnetosphere is not compressed. They obtained a magnetosphere radius where the 
part of the magnetic field of an NS outside rm is compressed according to the conservation of the 
magnetic flux at the accretion disk plane. The expression of the magnetosphere radius is 

= 5.54.cm , (23) 


where a is the viscosity parameter and ao.i = a/0.1, / = (1 


Aii/4 


Besides the above magnetosphere radius from theoretical modeling, Kulkarni & Romanova 
(2013) found from numerical simulation another dependence of the magnetosphere radius on the 
surface magnetic field, the mass, the radius of an NS and the luminosity of a XB. In their three- 
dimensional magnetohydrodynamic simulation of an accretion process at the quasi-equilibrium 
state, the gravitational, centrifugal, and pressure gradient forces were considered. It is different 
from Shi, Zhang &: Li (2014) in that the globular magnetic field is compressed by the accretion 
matter, i.e., the magnetosphere within the magnetosphere radius is also compressed. They obtained 
a magnetosphere radius as follows 


rm 3 = 7.12 * cm. 


^37 


"12 


(24) 


Please refer to Figure 2 in Shi, Zhang & Li (2014) for comparisons of the three models of magne¬ 
tosphere radius. 


2.2.3. The New Torque Model 

As seen from Section 2.2.1, the expression (Equation (20)) of the dimensionless torque will 
not be changed if we adopt the different expressions of the magnetosphere radius, but only a true 



- 10 - 


magnetosphere radius can lead to a correct dimensionless torque. The three kinds of magnetosphere 
radii in Subsection 2.2.2 are all admitted. In order to compare the magnetic field of an NS in BeXBs 
with the uncompressed magnetosphere to the one with the compressed magnetosphere below, we 
substitute the three expressions of magnetosphere radius (see Equations (22)-(24)) into Equation 
(20). Then we obtain 

rmi : (1 - w) + ^(| -2uj + w^) , 

(l)for uj<l, n = < ^-^2 : (1 - w) + 314.258 * / 34 /i 0 p-i/i 2 ^- 3 / 20 ^-i/i 2(2 _ 2oj + , (25) 

V m? . ■ (1 — cj) -|- 543.248 * ^^(l — 2aj + , 


rmi ; (1 - w) - 1 ) , 

( 2 )for a; > 1 , n = < rm 2 : (1 - w) + 314.258 * / 34 /i 0 p-Vi 2 ^- 3 / 20 ^-i/i 2 ( 2 ^-i _ ^ (26) 

_ rm 3 : (1 - w) -b 543.248 * Picj(|w-^ - 1) . 

As we can see in Eigure 2 , unlike the five-dimensionless torque models in subsection 2 . 1 , our 
result is convergent when the magnetosphere radius is close to the corotation radius (also see 
Equations (28) and (29)). Like some other dimensionless torque models, the zero torque acting on 
an NS is also not at the spin equilibrium state due to the assumption that the interaction between 
the magnetic field and the accretion disk also exists when n = rm = Vco, i.e., oj = 1 (see Equations 

(3), (23), and (24)). It should be noted that the dimensionless torque for r m 2 and rms do not 
depend only on the fastness parameter (w), but are also related to P (for r ^2 and rms) and L (for 
?’m 2 )) respectively. 

The dimensionless torques for the magnetosphere radii of the compressed magnetic fields 
(rm 2 & ?"m 3 ) change much more rapidly than the dimensionless torque for r^i and the reason 
can be found from Equation (27). The surface magnetic fields of NSs for rni 2 ^ ^’m 3 change with 

magnetosphere radius more rapidly than that for r^i (see Figure 2 in Shi, Zhang &: Li 2014). 

The values of the negative dimensionless torque from Equation (27) change only with when we 
compare the dimensionless torque in Figure 2 with the same fastness parameter. So the rapidly 
changed surface magnetic fields can lead to the large difference of the three panels of Figure 2. The 
last two panels of Figure 2 in fact express a relation between the surface magnetic field of NSs and 
the fastness parameters. 


3. Estimating the Surface Magnetic Fields of NSs in BeXBs 

Ghosh &: Lamb (1979) first calculated the magnetic moment of several XBs in the range 
3 * 1029 -4* 10^2 G • cm^ by their dimensionless torque model (GL). Klus et al. (2014) computed 
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Fig. 2.— Relation between the dimensionless torque and the fastness parameter for rmi, ?’m 2 ) and 
rjn 3 , respectively. 

the magnetic fields of NSs in BeXBs in SMC by GL and KR with the Alfven radius (i.e. r^i) and 
they obtained two kinds of solutions (the upper or lower solutions). Ho et ah (2014) argued that 
the upper solutions reflect the real magnetic field of NSs in BeXBs. 

Here we we extend the calculation on the surface magnetic fields of the NSs in BeXBs for 
Wang’s three-dimensionless torque models (Wang 1995) and the model (KR) in disk accretion. In 
order to compare the characteristics of the sources in disk accretion with that of the sources whose 
accretion modes are unknown, we calculate the the magnetic fields for all sources from Klus et ah 
(2014) except SXP 701 (H = 0 for SXP 701) with a disk accretion condition. Klus et al. (2014) 
computed the results for GL and KR with ta, but the inner radius was obtained as Equation (15), 
not TA in KR. Here we recalculate the surface magnetic held of NSs in BeXBs with the inner radius 
obtained by Equation (15). We adopt all the sources from Klus et al. (2014) except the source (SXP 
701) in subsections 3.1 and 3.2 and also the BeXBs in the Milk Way (MW), the Large Magellanic 
Cloud (LMC), and SMC in subsection 3.3 because we do not know the spin period derivatives of 
many sources in LMC and MW. 


3.1. Estimating Surface Magnetic Field With Four Available Torque Models 

In this work we select the characteristic values of NSs in accretion X-ray pulsars, i.e., M = 
1.4Mq, R = 10® cm, a = 0.1 (King et al. 2007) and so the inertial moment of a spherical NS is 
1.12 * 10^® g • cm^. From Equations (1), (3), and (5), we obtain 

= -2.72 * 10® * , (27) 

where the unit of P is s yr“^. Then the fastness parameter (tc) can be obtained when the dimen¬ 
sionless torque (Equation (13) or (14)) is substituted into Equation (27). 

We then compute the magnetic held solutions of NSs in BeXBs when Equations (7) and (10)- 
(12) are solved using = ^/r®, and the surface magnetic helds of NSs in BeXBs in Wang’s models 
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B= I 


(28) 


can be described as follows 

[ (A) : 6.05 * - w)-V2 Q, for w < 1 , 

6.05 * - W-1)-V2 G, for w > 1 , 

(B) : 6.05 * G , 

(C) : 1.56 * * 1%^ * ryiu^/^{l - w)-V3p-V6 ^ ^ 

where rco,8 = ?’co/(10® cm), pi =p/{l dyn/cm^). 

Similar to the above processes, we can obtain the surface magnetic field of NSs in BeXBs in 
KR with Equations (7), (14), (15), (22), and (27) 


B = 1.50 * 10^^LiyMiy^Rg * (i _ ^)-i/2 q 


(29) 


where Pi = P/1 s. 


As shown in Figure 3, we show the results of the NSs in BeXBs for Wang’s models and KR with 
their respective radius according to the different characteristic values of the fastness parameters. 
The detailed values of the fastness parameters for different sources are different, but the upper 
branch of the solutions in one model corresponds to a similar fastness parameter and we mark the 
mean values of the parameters in Figure 3. The fastness parameters for the lower branch of the 
solutions are very different and we only mark their scope (cj <C 1; the same process with Figure 
4). Most solutions with the fastness parameter (oj ~ 0.9) are higher than the critical magnetic 
field, indicating the same property of all BeXBs. Besides those solutions, there are several special 
solutions in all the panels for the source (SXP 4.78) with different fastness parameters and there 
is no magntic field solution for in GL (see Klus et al. 2014). There is no solution for the magnetic 
field for the sources (SXP 2.37, SXP 11.5, and SXP 16.6) in the above models. 


3.2. Estimating Surface Magnetic Field Using Our New Torque Model 

In this subsection we apply our new model presented in Section 2.2.3 and estimate the magnetic 
fields of NSs in BeXBs; note that here the magnetosphere radius outside the radius of an NS is 
considered as the inner radius of the accretion disk. 

Combing Equations (7) and (22)-(24), we can obtained the expressions of the fastness param¬ 
eters. After combing Equation (25)-(27) with the expressions of the fastness parameters, we can 
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1 10 100 1000 1 10 100 1000 


P (S) P (s) 

Fig. 3.— General results of Equations (28) and (29) between B and P for four-dimensionless torque 
models (Wang A, B, and C and KR). Hereafter in all figures R* indicates the radius of an NS, “up” 
indicates “spin-up” and “down” indicates “spin-down”; the dashed lines mark the quantum critical 
field, i.e., B = 4.414 * 10^^ G; the data for SXP 4.78 with a star. The filled symbols express the 
sources in disk accretion and the open symbols express the sources in unknown accretion modes, 
as identified in Klus et al. (2014). 


obtain the surface magnetic fields of NSs in BeXBs as follows 


Bu 


' 5.04 * 1014[(-P/n)-V2p,^ 

< r^2- 3.70* 10l5a();^/20^i7/io^4/5^203^40j^_p^^^_l/2p^^59/135j_27/4 ^ 

_ rm3 : 6.18* 10^2pf/4^11/4j^_p^^^_i/2p^^9/20j_io _ 


(30) 


As shown in Figure 4, our result shows the solutions of Equations (30) when the dimension¬ 
less torques (Equations (25) and (26)) are substituted. Similar to the results of the above hve- 
dimensionless torque models, the stronger magnetic held solutions in our model also have a similar 
fastness parameter. It means that the magnetic held has a relatively steady relation with the spin 
period when omitting the small change of the luminosity in BeXBs, suggesting that these solutions 
are near spin equilibrium. The fastness parameters for these stronger held solutions are lower than 
that obtained with the hve-dimensionless torque models in subsection 2.1, which means that the 
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equilibrium radii are smaller here. The lower solutions with different and very small fastness pa¬ 
rameters do not have a steady relation with the spin period, because these solutions are far from 
spin equilibrium. 

About the two branches of solutions in the left (for Tmi) and right panels (for Tma), the upper 
solutions are obtained for all these BeXBs, but the lower solutions are obtained only for some 
of them. There are two possibilities for this: (1) the upper and lower solutions give the correct 
magnetic fields for different BeXBs. However, in this case it is difficult to understand why the XBs 
have two very distinct groups of magnetic fields and there are almost no BeXBs with magnetic 
fields in between. (2) The lower solutions are non-physical and thus should be abandoned, and 
thus only the upper solutions are reasonable, as argued by Ho et al. (2014). In addition, all the 
lower solutions are for the spin-up NSs and have very small fastness parameters, i.e., far from the 
spin equilibrium. It was observed that some spin-up NSs in XBs sometimes also transit spin-down 
without significant luminosity change, implying they are not far from spin equilibrium. 

In fact, in the middle panel (for r^ 2 ) of Figure 4, only the upper solutions are obtained in near 
spin equilibrium, which indicates that the model of r ^2 is more reasonable than either or r^^,- 
We therefore conclude that the upper solutions with the r ^2 model give the correct magnetic fields 
of these NSs, which are also in near spin equilibrium. 

Note that although our new model presented above is an improvement beyond Klus el at. 
(2014), our conclusions are still affected by the simplistic assumptions about the magnetosphere 
radius calculations. Better magnetosphere radius models are still needed for realistically estimating 
more the magnetic fields when the detailed interactions between the disk and the magnetosphere 
are fully considered and the compression effect of the magnetic field in the accretion process is 
better described. 



Fig. 4.— Relation between the surface magnetic field of the NS (H 12 ) and PL^ for rmi, and 
?’m 3 , respectively. The filled symbols express the sources in disk accretion and the open symbols 
express the sources in unknown accretion modes, as identified in Klus et al. (2014). 
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3.3. Estimating the Surface Magnetic Field by Spin Equilibrium 


A spin-evolution timescale {P/P ~ 10^ — 10^ yr) of XB pulsars in the SMC (SXPs) can be 
obtained from the data of Klus et al. (2014), which is much shorter than the lO'^ year age of an 
OBe star and the decay timescale 10® yr of the global magnetic field in BeXBs (Ho et al. 2014). 
Therefore, BeXBs with strong magnetic fields in the SMC have enough time to achieve near spin 
equilibrium and the surface magnetic fields of NSs in most BeXBs will not be too low when they 
achieve the spin equilibrium. 


The above models are all including non-zero derivatives of the spin period {P ^ 0) and are 
called non-spin equilibrium models (Klus et al. 2014). We also can obtained the surface magnetic 
field of an NS with the spin equilibrium condition (rm = Vco) for rmi, r ^2 and as follows 


Bi2 


_ r^3 : 


(31) 


In order to compare the solutions from spin equilibrium with those from non-spin equilibrium, 
we show all the solutions (the solutions from Section 3.1 by non-spin equilibrium models, the 
solutions from the new model in Section 3.2 for r ^2 and the solutions from Equation 31 by spin 
equilibrium for r^ 2 ) in Figure 5. According to the criterion of disk accretion by Klus et al. (2014), 
we showed the solutions for the sources in disk accretion (urei/fcrei < 1 corresponding to the sources 
in disk accretion) in the left panel of Figure 5, and the solutions for the sources whose accretion 
mode has not been determined due to no available Urei in the right panel, where Urei is the relative 
velocity of accreted matter and Ucrei is its maximum relative velocity for which disk accretion can 
take place. We obtain the surface magnetic field of the NSs in the BeXBs in unknown accretion 
mode by the models of disk accretion. The result for the sources in wind accretion can be referred 
to Klus et al. (2014). As shown in Figure 5, the solutions with larger characteristic values for a; 
from the non-spin equilibrium condition are very close to the solutions from the spin equilibrium 
condition. 


Figure 6 shows the distribution of the spin frequency derivatives {v) of NSs in BeXBs in the 
SMC, which shows that the values of most v are close to zero and thus most BeXBs are quite close 
to spin equilibrium. A lot of observations on the sway of the spin period of an NS between spin-up 
and spin-down (e.g. Chakrabarty et al. 1995; Bildsten et al. 1997; Nelson et al. 1997; Postnov et 
al. 2015) suggest that the accretion processes of NSs in BeXBs are very close to spin equilibrium. 
It is consistent with the conclusion from the above analysis. 
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Fig. 5.— Composed figure of the surface magnetic fields of NSs from spin equilibrium and non¬ 
spin equilibrium conditions, in which the data come from Figure 3, our torque model for rjn 2 and 
from the disk spin eqnilibrinm condition for r^ 2 - Here “Eqnilibrinm” indicates the solntions from 
spin equilibrium, “1, 2” indicates the first and second solntions for some sources from non-spin 
equilibrium. The data marked with stars come from the special sources ( SXP 4.78). Left: the 
sources are in disk accretion. Right: the sonrces are in unknown accretion modes. 

Because most BeXBs are very close to spin equilibrium, we can compute the surface magnetic 
fields of the NSs in BeXBs in the MW, the LMC and the SMC in the spin equilibrium condition, as 
shown in Figure 7. We also show the surface magnetic field calculated from the observed cyclotron 
lines of some sources with B 12 = 0.863 * Eiq kevj if we suppose that the electrons are in the ground 
state and the line comes from the surface of an NS, where Eiq ^ev is the energy of the cyclotron 
absorption line in units of 10 keV. The result suggests 10^^ — 10^^ G surface magnetic fields 
of NSs in BeXBs, consistent to that reported by Caballero & Wilms (2012). If we consider the 
cyclotron line energy of protons in the same way, the surface magnetic field of an NS should be 
Bi 2 = 1586 * Eiq keV; which is much closer to the upper solutions discnssed above. Maybe the 
discussed cyclotron lines are indeed produced by protons. The result has also been discussed in the 
isolated magnetars by Harding & Lai (2006). 

The solutions of r ^2 and have much higher magnetic fields than those of r^i because of the 
compression effect of the magnetic field. If the compression effect of the magnetic field exists, the 
magnetic fields of most NSs in BeXBs will exceed the quantum critical magnetic field (44.14 TG), 
but all of them will not exceed the maximum “viral” value 1.77 * 10 ^®Mi. 4 Qi?^^ G, which is defined 
according to the condition that the magnetic energy ( ^^^^ * ■^) eqnals the gravitational binding 
energy ( ^^ ). It seems that there is a ceiling effect of the snrface magnetic fields of NSs in BeXBs 
dne to the maximum “virial” value as shown in Figure 7. 










- 17 - 




-16 -12 -8 -4 0 4 8 12 16 100150200 


Dot(v) 


Fig. 6.— Distribution of the spin frequency derivative of the NSs in BeXBs in the SMC. 


4. Discussion 
4.1. Wind accretion 

Although many observations indicate that there might be an accretion disk in a BeXB (Porter 
&: Rivinius 2003, James 2010), a system containing an OBe star with a non-truncated circumstellar 
disk might also be in the wind accretion state (Shakura et al. 2012). Shakura et al. (2012) ascribed 
the phenomenon of torque reversals to wind accretion; however, torque reversals have also been 
observed in disk-fed pulsars. Nevertheless, it is reasonable to assume that an accretion X-ray pulsar 
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Fig. 7.— Relation between the surface magnetic fields of NSs in BeXBs in and the spin period 
of NSs. The stars mark the magnetic field inferred from the observed cyclotron line sources (data 
from http://www.sternwarte.uni-erlangen.de/wiki/doku.php?id=cyclo:start). The data on the spin 
period of NSs and luminosity are from Knigge et al. (2011), Reig (2012), Kins et al. (2014), Cheng 
et al. (2014), and references therein. The circle points indicate the solutions for rms, the triangle 
points for rjn 2 j the square points for r^i, the solid points indicate the solutions from the sources 
for the SMC, the half-filled points for the LMC, and the open points for the MW as well. 
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swinging frequently at the zero torque should be in wind accretion state. Klus et al. (2014) found 
that the ratio Uj-ei/ucrei can be obtained for nearly half of the the systems in their sample of BeXBs 
and all of these systems are in disk accretion; the other half may contain some systems in wind 
accretion (also see Figure 5). 

However all our calculations presented here are only valid for the disk-fed pulsars. The models 
of Wang (1995), Kluzniak & Rappaport (2007) and Dai & Li (2006) were all developed only for 
disk accretion. For pulsars not in spin equilibrium, our dimensionless torque model evolved from 
Dai & Li (2006) and Wang (1995) cannot be extended to wind accretion; the main reason is that 
the model of magnetic field topology used by us is derived from disk accretion (Wang 1995). As 
a matter of fact, Klus et al. (2014) also did not calculate the surface magnetic field of the NSs of 
BeXBs not in spin equilibrium for wind accretion. Therefore, further development is still required 
to calculate the magnetic field of the NSs in BeXBs not in spin equilibrium for wind accretion. 

For pulsars in spin equilibrium, we assumed that the corotation radius is equal to the magneto¬ 
sphere radius for a disk-fed pulsar. However, this assumption is clearly in valid for wind accretion. 
Klus et al. (2014) assumed that the torques of spin-up and spin-down cancel each other for wind 
accretion and calculated their surface magnetic fields this way. The magnetic fields derived by Klus 
et al. (2014) from two spin equilibrium models on wind accretion are lower by approximately one 
or two orders of magnitude than that for disk accretion in spin equilibrium obtained by them or 
by us using the compressed magnetosphere, respectively. 

Is it then possible that the NSs in BeXBs in wind accretion have systematically lower surface 
magnetic fields than those in disk accretion? We believe this is unlikely because the surface mag¬ 
netic fields of NSs in BeXBs should be their intrinsic properties at birth and thus not influenced 
significantly by their different accretion modes and processes at a relatively small age in comparison 
with the decay timescale of the magnetic field. Given the relative reliability and simplicity of disk 
accretion models, the super strong magnetic fields calculated with disk accretion models might 
reflect the common characteristic of all NSs in BeXBs if a significant fraction of the BeXBs are 
indeed in disk accretion. This is consistent with the fact that the identified disk-fed systems cover 
the same parameter ranges as other systems, as shown in Figures 3-5. 


4.2. Super Strong Magnetic Fields 

Most young NSs have magnetic fields with orders of magnitude 10^^ G, but much lower mag¬ 
netic fields with orders of magnitude 10® G are always reported in LMXBs (e.g. Shi, Zhang & Li 
2014). In contrast, much higher magnetic fields with orders of magnitude 10^^ G are reported in 
magnetars, including anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs; Woods & 
Thompson 2006, p. 547; Mereghetti 2008). The most effects of the magnetars are ascribed to the 
strong magnetic fields, and they are especially powered by the strong magnetic fields. If a young 
NS (with a strong magnetic field 10^^ ~ 10^^ G ) has an evolutionary connection with the LMXBs, 
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the damping timescale of the magnetic field is 10® year (Harding &: Lai 2006). Thus, it seems that 
the lower magnetic field solutions with very small fastness parameters cannot be achieved because 
the lower magnetic fields are much lower than the fields of LMXBs and the evolution timescale 
should be much longer. We found that the surface magnetic helds of most NSs in BeXBs are higher 
than 10^^ G. Reig (2012) discussed the observation of the BeXB (4U 2206+54) and they thought 
that BeXBs are accreting magnetars. The NSs with strong magnetic fields in BeXBs as accreting 
magnetars may supply a channel to evolve into isolated magnetars because the two kinds of NSs 
have very similar characteristics, such as long spin periods and high spin period derivatives. 

The inferred extremely strong magnetic fields of NSs in some BeXBs supply a key clue to 
studying some basic physical problems (see Harding & Lai 2006). (1) Photon splitting becomes 
dominant as an attenuation mechanism, but photon splitting cascades may be prevented due to 
Adler’s selection rules (Harding & Baring 1996; Harding et ah 1998). Thus, that radiation from low 
energy photons in BeXBs should be absent or weak (e.g. radio radiation). (2) The equation of state 
of the bulk interior of an NS may be affected by the magnetic field since the field approaching the 
“virial” value. The NSs in those BeXBs would have lower density and the area of the radiation region 
would become larger than the NSs with weak magnetic fields because the high magnetic energy 
comparable with the gravitational potential energy compensates the interaction among baryons. 
The general relativity effect (such as gravitational redshift and the space-time curvature) would 
be alleviated because of the smaller mass-radius ratio. (3) Thomson & Compton cross sections 
are strongly reduced for the photon energies much lower than the cyclotron energy if the magnetic 
field is very strong (Paczynski 1992), and the maximum luminosity of an accreting BeXBs may 
exceed the Eddington luminosity. The ultraluminous X-ray sources (NuSTAR J095551+6940.8) as 
a accreting magnetar in HMXBs (Bachetti et al. 2014; Ek§i et al. 2015) may be an example of the 
BeXBs. 

5. Summary 

In this study we reviewed five-dimensionless torque models of X-ray binaries and presented 
a new model of the dimensionless torque based on the calculation of Dai & Li (2006). Using 
these models, we calculated the surface magnetic fields of the NSs in BeXBs in SMC. Finally, we 
estimated the surface magnetic field according to the spin equilibrium, including the NSs in BeXBs 
in the LMC and the Milky Way. Our conclusions can be summarized as follows. 

1. The dimensionless torque in the hve previous dimensionless torque models is divergent when 
the magnetosphere radius is equal to the corotation radius. 

2. The dimensionless torque in our new model is convergent and depends on the NS spin period 
and luminosity of X-ray binaries besides the fastness parameters, when applied to the two 
compressed magnetosphere radius models. 
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3. Two branches of solutions are found with all previous torque models, i.e., one with larger 
fastness parameters very close to spin equilibrium and the other with very small fastness 
parameters corresponding to spin-up far from spin equilibrium. The upper branch solutions 
with higher magnetic fields are close to spin equilibrium and more reasonable than that in 
the lower branch. 

4. By applying our new dimensionless torque with the compressed magnetosphere radius model 
in Equation (23) (Shi, Zhang & Li 2014), only one branch of solutions is obtained, corre¬ 
sponding to the upper solutions obtained in other models. This combination of dimensionless 
torque and magnetosphere models is thus favored. 

5. By assuming spin equilibrium, the estimated surface magnetic fields of the NSs in BeXBs in 
the SMC, the LMC and the Milky Way show similar characteristics. 

6. The estimated surface magnetic fields for the compressed magnetosphere models are much 
higher than that for the uncompressed magnetosphere model. They exceed the quantum 
critical magnetic field, but do not exceed the maximum “virial” value. 
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